I 


NASA Contractor Report 181883 
ICASE Report No. 89-55 





f 


ICASE 


ON THE LARGE-EDDY SIMULATION OF 
TRANSITIONAL WALL-BOUNDED FLOWS 


Ugo Piomelli 
Thomas A. Zang 
Charles G. Speziale 
M. Y. Hussaini 


Contract No. NAS 1-1 8605 
August 1989 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

Operated by the Universities Space Research Association 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665-5225 

( NASA— CR— 1 81883) ON THE lARGE-POOY N39-2B763 

SIMULATION OF TRANSITIONAL WALL-BOUNOFO 

FLOWS Final Report (NASA. Lanqley Research 

Center) 29 p CSCL 200 Unci as 

HI/34 0231751 


ON THE LARGE-EDDY SIMULATION 
OF TRANSITIONAL WALL-BOUNDED FLOWS 


Ugo Piomelli* 

Department of Mechanical Engineering 
University of Maryland 
College Park, MD 20742 

Thomas A. Zang, 

NASA Langley Research Center 
Hampton, VA 23665 

Charles G. Speziale and M. Y. Hussaini* 
ICASE, NASA Langley Research Center 
Hampton, VA 23665 


Abstract 


The structure of the subgrid scale fields in plane channel flow has been studied at var- 
ious stages of the transition process to turbulence. The residual stress and subgrid scale 
dissipation calculated using velocity fields generated by direct numerical simulations of the 
Navier-Stokes equations are significantly different from their counterparts in turbulent flows. 
The subgrid scale dissipation changes sign over extended areas of the channel, indicating 
energy flow from the small scales to the large scales. This reversed energy cascade becomes 
less pronounced at the later stages of transition. Standard residual stress models of the 
Smagorinsky type are excessively dissipative. Rescaling the model constant improves the 
prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. 
Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results 
of a large eddy simulation of transition on a flat-plate boundary layer compare quite well 
with those of a direct simulation, and require only a small fraction of computational 
effort. The inclusion of non-dissipative models, which could lead to further improvements, 
is proposed. 


‘This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1. Introduction 


The transition from laminar to turbulent flow has been the subject of much research. 
Experiments, theoretical studies and numerical methods have been brought to bear on the 
subject. Direct numerical simulations (DNS) of the Navier-Stokes equations are particularly 
useful for the study of the late stages of transition, in which nonlinear effects become signif- 
icant and linear theory as well as perturbation methods fail. This technique has been used 
to study the structure of the detached shear layer and lambda vortices that arise due to the 
secondary instability, 1-4 and the effects of heating, 5 streamwise vortices 6 and curvature. 7 
Recently, Gilbert and Kleiser 8 have carried out the direct simulation of channel flow at a 
Reynolds number Re c = 5000 (based on centerline velocity and channel halfwidth) for the 
full transition process from laminar to turbulent flow. 

Direct simulation of the governing equations has so far been limited to simple flows and 
low Reynolds numbers due to the massive computational effort it requires; the use of two 
million grid points is not uncommon, and simulations have been performed with up to five 
million points. 9 In general, more resolution is required to study the later stages of transition 
than for a turbulent flow, since in the former case it is essential to distinguish between truly 
chaotic behavior and grid-scale fluctuations caused by insufficient resolution. This resolution 
requirement makes the application of direct simulation to the study of flows of engineering 
interest unfeasible, at least in the near future; less computationally intensive methods are 
required for engineering applications. 

A technique that has been successfully applied to the study of turbulent flows is large- 
eddy simulation (LES). In LES one accurately resolves only the large, energy-carrying scales 
of the motion and models the effect of the small scales (subgrid scales), which appears 
in a residual stress term. In turbulent flows, the small scales tend to be more isotropic 
and homogeneous than the large scales, and do not depend very strongly on the boundary 
conditions. For these reasons their effect can be represented by simpler models than the 
currently used turbulence models. Large-eddy simulations have been successfully applied to 
a variety of turbulent, wall-bounded flows such as plane channel flow, 10-12 boundary layers 13 
and channel flow with transpiration. 14 The only applications of this technique to transitional 
flows to date are due to Horiuti, 15 Dang and Deschamps 16 and Deschamps and Dang. 17 

Horiuti 15 used a simple residual stress model to extend his direct simulations to the later 
stages of transition. The model that he employed, the Smagorinsky 18 model, dissipates a 
significant amount of turbulent energy from the resolved scales, effectively slowing down 
the transition process. Horiuti’s results are, therefore, only qualitatively correct, due to the 
significant inaccuracies of the residual stress model. 

In two related papers Dang and Deschamps 16 and Deschamps and Dang 17 describe their 
large-eddy simulation of transitional flow. Their results are inconclusive: they do not make 
any comparison with theoretical predictions at the early stages of transition, and, at the later 
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stages, compare their LES results with those of a direct simulation that they also performed. 
Their DNS results, however, compare quite poorly with experimental data, perhaps due to 
insufficient resolution, and no conclusive statement can be made regarding the accuracy of 
their large-eddy simulation. 

It is the purpose of the present paper to establish a rational basis for the application of 
LES to transitional flows. This will be done by examining the results of direct simulations of 
transitional channel flow to understand the physical behavior of the subgrid scale structures. 
The energy exchange between subgrid scales and large scales will be studied and compared 
with turbulent flow behavior. It will be shown that the interaction between large and small 
scales in transitional flows is quite different than in turbulent flows. 

In Section 2 the governing equations will be laid out and basic considerations for the 
optimal choice of filter for the LES of transitional flows will be made. In Section 3, the 
direct simulation databases will be examined in detail, and residual stresses and subgrid 
scale dissipation will be calculated and presented. Commonly used residual stress mod- 
els will be evaluated in Section 4, and preliminary results of a large-eddy simulation of a 
transitional flat-plate boundary layer flow will be discussed in Section 5. Conclusions and 
recommendations for future work will be made in the last section. 


2. Numerical formulation 


In large-eddy simulations we decompose the flow variables (velocity U{ and pressure p ) 
into a large scale component (denoted by an overbar) and a subgrid scale component. The 
large scale field is defined by the filtering operation: 


f{^x\ } X2,X3^ — / £?i(x,', a , x 2) x^Jdx^dx2clx^ 

JD *=1 


( 1 ) 


where the integral is extended over the entire domain D and Gi is the filter function in the 
zth direction. 19 

The filtered Navier-Stokes and continuity equations, which describe the evolution of the 
large, energy-carrying eddies, can be obtained by applying the filtering operation to the 
incompressible Navier-Stokes and continuity equations to yield 


dui 

dt 



dui 

dxi 


dp dTij 1 d 2 u l 

dxi dxj Re dxjdxj 

0, 


( 2 ) 

( 3 ) 


in which a reference length and velocity scale are used to make U{, p, Xi and t dimensionless. 
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The effect of the subgrid scales appears in the residual (or subgrid scale) stress 


Tij = U t Uj - UiUj. 


( 4 ) 


The residual stress can be further decomposed into the cross stress Cij and the subgrid scale 
Reynolds stress Rij, given by 


Cij = UiU'j + UjU'i 

Rij = <u'j 


(5) 

( 6 ) 


respectively, in which u[ = rq — U{ is the subgrid scale velocity. In addition, the Leonard 
stress Lij is defined as 

Lij = W&j - U iUj. (7) 

The cross stress Cij and subgrid scale Reynolds stress Rij are invariably modeled, while the 
Leonard stress Lij can be either computed or modeled. 

Most of the models currently in use correlate the residual stress with the large scale 
strain-rate tensor Sij, 

1 / f)Tij \ 

( 8 ) 


- _ 1 fduj duj 

13 2 \dxj dx ij 


which, in incompressible flow, is traceless. The trace of the residual stress tensor is, therefore, 
commonly incorporated into the total pressure, and only the anisotropic part of the residual 
stress tensor 

T ij = T H ~ S H T kk/3, (9) 

(in which r = q 2 sgt is the subgrid scale kinetic energy and Sij is Kronecker delta) requires 
modeling. 

For the large-eddy simulation of turbulent flows, the most commonly used filters are 
the Gaussian filter, the sharp Fourier cutoff filter and the box filter. The Gaussian filter is 
defined as 


( 6 \ 1/2 

Gi(s,-,x') = f exp 


-6(xi - x') 2 


A? 


( 10 ) 


where A,- is the filter width in the ith direction, usually related to the grid size Ax^. 
The sharp Fourier cutoff filter is most conveniently defined in Fourier space: 


Gi(ki) = 


1 for ki < Kd 
0 otherwise 


( 11 ) 


in which Gi(ki) is the Fourier coefficient of the filter function Gi, and Kd is the cutoff 
wavenumber in the zth direction, that is related to the filter width by Kd = 7t/A 
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The variable-width box filter is defined as 


G;(Xj, X ( ) 


1/A; for |x; — x(| < Ai/2 
0 otherwise. 


( 12 ) 


All of the filters defined above have been used for the large-eddy simulation of turbulent 
wall-bounded flows. 

When either the Gaussian or the box filter is used, the subgrid scales must represent 
both large and small scales of motion. 3 * * * * * * * * 12 * * * * * When the cutoff filter is used, on the other hand, 
the subgrid scales only contain the effect of the structures with wavenumber larger than 
the cutoff wavenumber K This characteristic of the cutoff filter is very desirable for the 
large-eddy simulation of transitional flows, in which the large scales of motion (such as the 
Tollmien-Schlichting waves) must be resolved very accurately, and large eddy simulation 
should give the same results as direct simulation in the early stages of transition. 

It is also useful to remark that the sum of the cross stress, Leonard stress and subgrid scale 
Reynolds stress is invariant to a Galilean transformation, regardless of the filter used. 20 The 
subgrid scale Reynolds stress Rij has the same property. The Galilean invariance properties 
of Leonard and cross stresses, however, depend on the type of filter used. If either the 
Gaussian or the variable-width box filter is employed, neither the Leonard stresses nor the 
cross stresses are Galilean invariant, but their sum is. On the other hand, if the cutoff filter 
is used, both Leonard and cross stresses are Galilean invariant. 


3. Analysis of databases 

To gain a better understanding of the structure of the residual stress tensor in transitional 

flows, one can use the results of direct simulations to evaluate the residual stresses, and study 

their behavior at various stages of transition, by comparing the exact residual stress, obtained 

from (9) and (4), with the predictions of a model. This technique, which is known as an a 

priori test, was first applied by Clark, Ferziger and Reynolds 21 to the study of homogeneous, 

isotropic turbulence. To calculate the exact residual stress, the velocity field obtained from 

a direct simulation of the Navier-Stokes equations is filtered using one of the filters defined 

above to yield the large scale velocity field rT t -; the subgrid scale velocity field u( can next 

be obtained by subtracting tT; from u;. The large scale and subgrid scale velocities can then 

be used to calculate the exact values of quantities of interest such as the residual stress r,y 

[which can be obtained from (4)] and the dissipation by the residual stresses 

= TijJij = r^Sij. (13) 

This term, which appears as a dissipation term in the transport equation for the resolved 

energy and as a production term in the transport equation for the subgrid scale energy, 
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represents the energy transferred from large scales to subgrid scales (the “energy cascade” 
in turbulent flow). It has been argued 22 that, as long as a subgrid scale model drains energy 
from the large scales at the correct rate, its overall performance will be acceptable even 
if the model fails to predict the local stress accurately; the considerable success obtained 
by large-eddy simulations which use the Smagorinsky 18 model has been attributed to this 
characteristic. Finally, the filtered field rq can be used to compute the predictions of the 
various models, which can then be compared with the exact quantities. 

Although the a priori test gives useful indications regarding the behavior of subgrid-scale 
quantities, its results should not be considered final. The ultimate test of the accuracy of a 
residual stress model is a posteriori, when the model is used in a large-eddy simulation and 
its results are compared with experimental and numerical results. 

In the present work we have used numerical databases similar to those generated by Zang 
and Krist 9 to study the behavior of the components of the residual stress tensor and of the 
subgrid scale dissipation during transition in incompressible plane channel flow (a schematic 
representation of the channel geometry is shown in Figure 1). The velocity field from the 
direct simulation, which had been obtained at Re c = 8000 was filtered using the cutoff filter 
in the xz-plane; no filtering was applied in the normal direction. The simulation was started 
with a parabolic velocity profile, on which a 2D and 3D disturbance were superimposed. 
The 2D disturbance was a Tollmien-Schlichting wave with a wavenumber a x = 1 (made 
dimensionless by the channel halfwidth 5) and an amplitude of 2%. The wavenumbers of 
the 3D disturbance were respectively a x = 1 and a z = 1 . 5 ; its amplitude was 0.02%, and 
its phase was matched with that of the 2D wave. The mass flux through the channel was 
maintained constant and equal to its laminar value. Channel flow databases where preferred 
to boundary layer databases for the a priori test because channel flow can be computed far 
more economically than boundary layer transition. 

As a preliminary step, the effect of the filter width on the a priori test results was 
studied. It was found that changing the filter width has two principal effects: in the first 
place, the amount of energy contained in the subgrid scales changes, as does their percentage 
contribution to the residual stress tensor. Moreover, removal of particular waves may result in 
slightly different residual stress profiles. The qualitative behavior of subgrid scale statistics, 
however, is not much influenced by the filter width; for this reason in the remainder of this 
study we will only present results obtained with a single filter width, corresponding to cutoff 
wavenumbers K c \ = K c 3 = 9. 

Five flow realizations were examined at times t = 175, 180, 185, 190 and 196 (made 
dimensionless by the centerline velocity U c and channel halfwidth). This covers the period 
from the development of the so-called “second spike” (at t — 175), through the start of 
the laminar breakdown (at t = 185), into the transitional region ( t = 196). The transition 
process itself is not complete until t — 240. The Reynolds stress < u"v" > (where u " = Ui~Ui 
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is the turbulent fluctuation in the fth direction, U{ =< tq >, and < . > represents an average 
over the xz- plane) reaches a maximum at t = 210. By t — 196, < u"v" > is approximately 
two thirds of its maximum value, and is nearly twice the final turbulent value. In all cases 
the velocity, which had been computed on grids using up to 216x192x216 grid points, was 
filtered onto an 18x192x18 grid. 

In Figure 2, the temporal histories of the subgrid scale kinetic energy q] g , are shown, 
along with its ratio to the total turbulent kinetic energy q 2 = u"u". At the late stages of 
transition the subgrid scale energy increases rapidly, as small scale structures are generated. 
For t > 180, subgrid scales account for almost 30% of the total energy. This result indicates 
that a finer mesh than the one used for the a priori test is required to resolve the late stages 
of transition. 

The (l,2)-component of the residual stress tensor is shown in Figure 3; the residual stress 
is less than 20% of the total stress (except in the near- wall region). For t < 180, the residual 
stress changes sign at approximately y/6 = —0.8; this behavior is not observed in turbulent 
flows. At later times the sign change present at the early stages disappears and, as small 
scale structures are generated in the near-wall region, the residual stress increases there. 

The sign reversal in the residual stress suggests that the energy cascade may be reversed, 
at least locally, in transitional flows: energy may flow from small scales to large scales. 
This hypothesis is supported by an examination of the subgrid scale dissipation (Figure 4). 
Although the mean dissipation is always negative (indicating that, on the average, the small 
scales absorb energy from the large ones) the fact that, for t < 180, the root-mean-square 
dissipation is approximately five times as large as the mean indicates that the reversed energy 
flow is significant in considerable regions of the channel. To clarify this point further, the 
the area of each xz-plane over which e, gs is positive, A + was calculated, and its ratio to the 
total area of an xz-plane, A t , is shown in Figure 5 at t = 175. In almost half of the channel, 
at this time, energy is flowing from the small scales to the large ones. At its peak, the 
rms dissipation by the subgrid scales is approximately 65% of the mean viscous dissipation 
e v = —2 SijSij/ Re, and 60% of the rms viscous dissipation. As transition progresses, however, 
the subgrid scales become more dissipative, especially in the near-wall region. At t = 196 
the rms c sgs is only 80% larger than the mean, indicating that the regions of reversed energy 
flow are much reduced in extent and intensity. 


4. Modeling considerations 

The most commonly used residual stress models for large-eddy simulations of turbulent 
flows are modifications of the Smagorinsky 18 model for the anisotropic part of the residual 
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stress. The Smagorinsky model is an eddy viscosity model which relates the modeled residual 
stress Tfj to the strain-rate tensor Sij through an eddy viscosity vt : 

t% = -2«/ T £; (14) 

in which v ? is given by 

v T = e^iSaSij ( 15 ) 

and l is a subgrid length scale, that is usually related to the filter width by 

l — Cs ( AxAyAz) 1 ^ 3 , (16) 


where the Smagorinsky constant Cs is set equal to 0.1 for wall bounded flows. 

The Smagorinsky model is absolutely dissipative. Computed subgrid scale dissipations 
are shown in Figure 6, and their integrated values over the entire computational domain 
are given in Table I. The model is incapable of predicting the reverse cascade regions, and 
overestimates the subgrid scale dissipation. Part of the problem is due to the fact that the 
model is only based on large scale quantities, and predicts nonzero residual stresses even in 
laminar flow. This limitation can be overcome in many ways. One can, for example, scale 
the modeled stress according to the stage of transition the flow is undergoing; the constant 
C s , for example, can be multiplied by the factor (Hi — H)/(Hi — H t ), where H — 8* / 9 is the 
shape factor, 8* is the displacement thickness, 9 the momentum thickness, and the subscripts 
l and t refer respectively to laminar and fully developed turbulent flow. For laminar flow, 
Hi — 5/2, while, for Re c = 8000, H t ~ 1.7. 23 The length scale l can then be given by 

t = Cs^~(AxAyAz)' 13 . (17) 

In turbulent flows the length scale is often multiplied by a van Driest damping function 24 
to account for the fact that the growth of even the large structures is inhibited by the 
presence of the wall. The length scale l then becomes 


i = C s 


1 — exp(— y + /25) (AxAyAz) 1 ^ 3 


(18) 


where y + is the distance from the wall made dimensionless by the shear velocity u T = 
( t w /pY I 2 (where t w is the wall stress) and kinematic viscosity v. Although this approach 
yields non-zero residual stresses even in laminar flow, these stresses are very small due to 
the small value of the shear velocity (and, therefore, of y + ) in laminar flow. In laminar flow, 
the subgrid scale dissipation is 35% of the viscous dissipation if (16) is used, zero if (17) is 
used and 7% of the viscous dissipation if (18) is used. 

Although the introduction of damping improves somewhat the prediction of the inte- 
grated dissipation (see Table I), the reversed energy flow is still not permitted, and the 
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dissipation is not well predicted locally. At later stages, all models underpredict the subgrid 
scale effects. The use of damping would probably lead to a more realistic prediction of the 
early transition process than the straightforward Smagorinsky model used, for example, by 
Horiuti. 15 However, the time development of single structures would probably still be incor- 
rectly predicted, since their growth rate would be damped where the model is excessively 
dissipative, and amplified where the model is not dissipative enough. 

Two models that have the capability of generating positive dissipation (i.e., reversed 
energy flow) are the mixed model 25 and the nonlinear eddy viscosity model. 26 The mixed 
model is commonly written as 


Tij = C B ( u { Uj - UiUj) - 2u T Sij (19) 

and has been successfully applied to the large-eddy simulation of turbulent flows. 12,25 The 
model in this form cannot be used with the cutoff filter, which has the property that U{ = f?;. 
A modification of this model, which can be used with the cutoff filter, is 


T- G B (^i ^i) (Uj tij) 2 V'B S 


(20) 


in which represents the velocity filtered using a sharp cutoff filter with a cutoff wavenumber 
K' ci < K c i ; here the Galilean invariance constraint is satisfied for any value of C E - 
The nonlinear eddy viscosity model 26 takes the form 

= - 2 q&ij - c D e (s ik s kj - ^s mn s mn s^ - c E e (f {j - ( 21 ) 

o 

where q is a velocity scale and Sij is the Oldroyd derivative of S{j‘. 



dS ij _ dSij 

r~ + Uk- 


dt 


dxi 


dui ^ duj ^ 
8x t ili dx k Sli 


(22) 


This nonlinear model has never been applied to large-eddy simulations of turbulent flow, 
although it gave more accurate predictions than the commonly used models when employed 
as a standard turbulence model for wall-bounded flows. 26 

Both the mixed model and the nonlinear model are not strictly dissipative, so that their 
use could lead to improved results. They, however, require added empirical information 
in the form of additional constants and the new cutoff wavenumber K which must be 
optimized for turbulent and transitional flows. 

The a ‘priori tests do not allow us to estimate to what extent the excessive dissipation 
affects growth rates and the final results, or whether the inclusion of energy-producing terms 
in the residual stress model is necessary. For this reason, large-eddy simulations which employ 
various variants of the Smagorinsky model have been carried out to evaluate a posteriori the 
accuracy of the model. These simulations will be discussed in the next section. 
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5. Large-eddy simulation results 


To evaluate more conclusively the effect of dissipation by the subgrid scales on the tran- 
sition process, two large-eddy simulations have been carried out in which the Smagorinsky 18 
model has been used to parameterize the residual stresses. In these simulations the transi- 
tion process on a flat-plate boundary layer was calculated. The flat-plate boundary layer was 
preferred over plane channel flow because of the greater similarity this flow has to problems 
of technological importance in aerodynamics. 

The simulations were performed using a spectral collocation method. Fourier expansions 
were used in the streamwise and spanwise directions, Chebyshev polynomials in the wall- 
normal direction. The code uses the improved splitting method described in Zang and 
Hussaini. 27 The time advancement is performed using an implicit Crank-Nicolson method 
for the diffusion term, residual stresses and incompressibility constraint; a third order Runge- 
Kutta method is applied to the remaining terms. 

The number of grid points used in this simulation was 24x48x24. Initial conditions 
for the simulation consisted of a Blasius velocity profile, on which a 2D and 3D disturbance 
were superimposed. The 2D disturbance was a Tollmien-Schlichting wave with a wavenumber 
a x = 0.25 (made dimensionless by the displacement thickness 6*) and an amplitude of 1.8%. 
The wavenumbers of the 3D disturbance were respectively a x = 0.25 and a z = 0.209; its 
amplitude was 0.8%, and its phase was 36° ahead of the 2D wave. The thickness of the 
boundary layer was assumed to be constant in the streamwise direction. The filtered Navier- 
Stokes equations (2)-(3) were integrated for a total time t/t^s = 4.5 (made dimensionless 
by the period of the Tollmien-Schlichting waves trs)- Large-eddy simulation results were 
compared with those of a direct simulation which used up to 144x144x324 grid points. 28 
The direct and large eddy simulation codes were essentially the same; the CPU time required 
per iteration and grid point was, therefore, the same. Due to the coarser grid used by the 
LES, however, less than 4 CPU hours were required to advance the solution for 4.5 Tollmien- 
Schlichting periods. The direct simulation, by comparison, required over 100 CPU hours to 
integrate the equations for the same period. 

The first simulation used the Smagorinsky model (14) with the length scale (18), based 
on van Driest damping. The dissipation by the subgrid scales was found, in this case, to 
be excessive. The growth rate o> obtained from the large-eddy simulation is compared in 
Figure 7 with that of the direct simulation. The growth of the perturbation is significantly 
damped in the LES, indicating that the effect of the residual stress model is too large. 

A new length scale combining Eqs. (17) and (18), was then used: 


i = C s ~~§-(AxAyAz) 1/3 

hi — h t 


1 - exp(— y + /25)] 


(23) 
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in which Hi = 2.6 for the Blasius profile, and H t was set to 1.4 on the basis of experimental 
evidence. 29 When this length scale is used, the model is essentially turned off during the 
early stages of transition; this allows for much more accurate prediction of the early stages 
of transition. The growth rate predicted by an LES using (23) is compared in Figure 7 with 
the DNS results. The agreement between the DNS and the LES is good at the early stages. 
For 2 < t/txs 5: 3.5 the model appears to dissipate excessively, slowing down the evolution 
of the perturbation. At later times, as the subgrid scales become more dissipative, the LES 
predicts an increased growth rate similar to that observed in the DNS results. This increase, 
however, occurs at a later time, perhaps due to the errors at the intermediate stages of the 
transition process. 

The development of the shape factor H (see Figure 8) follows similar trends: if no scaling 
is used, H remains at its laminar value throughout the simulation. When scaling is applied 
in the early stages, the LES results compare well with the DNS ones. The decrease in H that 
takes place when the mean velocity profile becomes significantly altered by the perturbation 
velocity occurs at a later time in the LES than observed in the DNS results, possibly because 
of the poor prediction of the intermediate transition processes. 

The harmonic energy content of selected scales is shown in Figure 9. The time evolution 
of the larger scales of motion (the 2D Tollmien-Schlichting waves are the (1,0) mode, for 
example) is predicted with surprising accuracy, although, once again, the transition process 
is delayed by the excessive dissipation. The large scale Reynolds stresses < u"u" > (where 
u" = Ui — Ui is the fluctuating large scale velocity) are also predicted with acceptable accuracy 
(see Figure 10). It should be noticed that, at this stage and with the mesh used in the present 
simulation, the subgrid scales account for about 20% of the total turbulent kinetic energy. If 
one were to compare the LES results with experimental data or with unfiltered DNS results, 
one should expect that discrepancy. In the present case, the DNS data was filtered onto a 
24x48x24 grid to yield a more significant comparison. 

Equal du/dy contours are shown in Figure 11. Although the large-eddy simulation cap- 
tures all the essential features of the flow, the liftup of the detached shear layer is not as 
pronounced, consistent with the delay of the transition process due to the excessively dissi- 
pative nature of the residual stress model. 


6. Conclusions 

A detailed study of the structure of the subgrid scales in transitional plane channel flow 
has been conducted at various stages of the transition process. It has been demonstrated that 
the residual stress and subgrid scale dissipation obtained from direct numerical simulation of 
the Navier-Stokes equations are significantly different from their counterparts in turbulent 
flows. During transition the subgrid scale dissipation changes sign for considerable areas 
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of the flow, indicating a reversal of the energy cascade wherein large scales absorb energy 
from the small ones. It was found that the standard residual stress models are excessively 
dissipative; they erroneously delay the onset of transition in the large-eddy simulation of 
boundary layers, leading to considerable error. The use of a shape factor function combined 
with van Driest damping to form a modified Smagorinsky model for transitional flows was 
shown to partially alleviate this problem. This new residual stress model was demonstrated 
to yield results for the disturbance growth rate and energy content of selected harmonics 
in transitional boundary layers that compared favorably with the direct simulations. These 
LES results were obtained at a small fraction of the computational expense required for the 
direct simulation. 

The simulations described above were carried out with the intent of determining the 
feasibility of large-eddy simulations for transitional flows. The fact that the LES results 
compared reasonably well with DNS results, indicates that this technique can be successfully 
applied to the simulation of transitional as well as turbulent flows. Since these results were 
obtained using very simple residual stress models, it appears that in transitional flows (as 
well as in turbulent flows) accurate calculation of the large scales is of primary importance 
for the prediction of quantities of engineering interest. 

While the results of this study have established the feasibility of large-eddy simulations 
with relatively simple residual stress models, there is still a need to develop improved models. 
For example, it may be useful to introduce an intermittency factor based on some large scale 
kinematic quantity that senses the development of small scale turbulence activity. In this 
manner, the residual stress model can be turned on or off depending on the local development 
of turbulence - a feature that could considerably improve the accuracy of transitional LES. 
The inclusion of non-dissipative models may also be of assistance. Two such models have 
been proposed: one based on a nonlinear eddy viscosity model and another based on the 
commonly used mixed model. The calibration and testing of such new models is a substantial 
research effort that is beyond the scope of the present paper. These, as well as other technical 
issues, are part of an ongoing research effort on the large-eddy simulation of transitional flows. 
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Figure 2. Subgrid scale energy and energy ratio obtained from direct simulations; o t = 
175; + t = 185; o t = 196. (a) q ] ga ; (b) q%,/q 2 . 
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Figure 6. Subgrid scale dissipation e sgs ; o DNS; A modeled, length scale given by 
(16); + modeled, length scale given by (17); x modeled, length scale given by (18). 
(a) t = 175; (b) t = 185; (c) t = 196. 
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Figure 8. Time evolution of the shape factor H. LES, length scale given by (18) 

LES, length scale given by (23); □ DNS. 
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Figure 9. Time evolution of the energy content of selected harmonics. The length scale 

for the LES was given by (23). (1,0) mode: LES, o DNS; (1,1) mode: LES, 

A DNS; (2,0) mode: LES, + DNS; (2,1) mode: LES, x DNS; (2,2) mode: 

LES, o DNS. 
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Figure 10. Large scale Reynolds stresses < u"u" >. The length scale for the LES was 
given by (23). The LES results are at tjt TS = 4; the DNS results are at t/t TS = 3.88 
and are filtered on a 24x48x24 grid, u" : LES, o DNS; v : LES, + DNS, 

w" 2 : LES, A DNS; u"v": LES, x DNS. 
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